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1 Introduction 



Some of the most precise determinations of the weak mixing angle sin 2 9 e ff come from 
measurements of asymmetries in fermion production on the Z peak [|IJ. In particular, the 
forward-backward asymmetry of b quarks is measured with a precision of about 2%, allow- 
ing an extraction of sin 2 9 e ff with almost per mille accuracy. However, since we are dealing 
with quarks in the final state, we must ensure that QCD corrections, both perturbative and 
non-perturbative, are understood to at least the same precision. From simple power count- 
ing, it is clear that this necessitates including 0(a\) perturbative and 1/Q non-perturbative 
effects. Even these will probably not be enough in the future, when linear e + e~ colliders 
are hoped to reach a precision of order 0.1% @. 

The 0(as) perturbative corrections were first calculated in Ref. || in the massless 
approximation. The mass corrections to this result were first calculated in Refs. 0] and 
were found to be significant ~ asmb/Mz- These calculations used a slightly different 
definition of the asymmetry than the experimental measurements, which use the thrust 
axis rather than the quark direction. This difference was rectified in Refs. || . 

To date there have been two 0(a\) calculations, both in the massless approximation 
using the quark direction. The classic calculation of Altarelli and Lampe J7| determined the 
0(«g) coefficient numerically and found it to be small. This result has been the basis of all 
the experimental analyses since. However, the recent analytical calculation by Ravindran 
and van Neerven || obtained a coefficient about four times bigger. This discrepancy is 
comparable to the size of the experimental errors and needs to be resolved before the 
final electroweak fits to the LEP1 data can be made. The C(ag)-calculation using the 
experimentally-used thrust axis definition, would also be highly desirable. 

In this paper we perform a numerical calculation of the 0(a\) corrections to the 
forward-backward asymmetry, and compare our results with the existing calculations. We 
also calculate for the first time the corrections using the thrust axis definition rather than 
the quark direction. 

The paper is set out as follows. In Sect. ^| we define the forward-backward asymmetry 
and the closely-related left-right forward-backward asymmetry @ and recall some features 
of the tree-level and 0(as) perturbative calculations. In Sect. [5] we discuss the general set- 
up of the C(a|) calculation, and divide it into several parts. We pay particular attention 
to the four-6 final state, which will turn out to play an important role in our calculation. 
In Sect. f| we make some final remarks on the details of the calculation, before presenting 
our results for the C(a|) coefficients with the two axis definitions. We also compare our 
results with the existing calculations. We discuss the impact of our results in Sect. |!| and 
try to estimate the remaining theoretical errors. We leave some more technical details of 
the calculation to Appendices A and B. 
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2 Definition and perturbative calculation 



The simplest definition of the 6-quarkQ forward-backward asymmetry Afb is 

_ N F -N B 

where Np and Nb are the number of b quarks observed in the forward and backward 
hemispheres, respectively. 

The axis that identifies the forward direction can be defined in a variety of ways. How- 
ever, for the purpose of making Afb computable in QCD perturbation theory, the axis 
must be defined in an infrared- and collinear-safe manner. In this paper we explicitly con- 
sider two different definitions: the 6-quark direction, and the thrust axis direction. The 
thrust axis has a two-fold ambiguity: we use the one that is nearer the 6-quark direction. 
In the following, the forward-backward asymmetries with respect to the 6-quark direction 
and to the thrust axis direction are denoted by A b FB and A FB , respectively. 

According to the definition in Eq. (|l|), A F b can be expressed in an equivalent way in 
terms of the cross section 

cKe+e- -> b + X ) (2 ) 
dx dcos 9 

for inclusive 6-quark production, where x is the fraction of the electron energy carried by 
the b quark and 6 is the angle between the electron momentum and the direction defining 
the forward hemisphere (both energies and angles are defined in the centre-of-mass frame). 

Starting from the distribution in Eq. fl2|), we can introduce the forward and backward 
cross sections op and ag: 

f 1 f 1 dcr f° f 1 dcr 

a F = dcos6» / dx- — -, a B = / dcos6» / dx- — -, (3) 

Jo Jo dxdcosfc 1 J-i Jo dzdcosfl 

and the symmetric and antisymmetric cross sections as and &a- 

(?S = (TF + <yB, Oa = (?F-<JB- (4) 

We can then write the forward-backward asymmetry as 

A F b = ~~~ ■ (5) 

o"5 



In the perturbative QCD calculation of as and a a, we have to evaluate the correspond- 
ing matrix element squared, which is given by the product L^ V T^ U of the leptonic and 
hadronic tensors L^ v and T^ v . Then we could perform the integration over the final-state 
parton momenta in T^ u and finally the integration over the scattering angle 9. Nonethe- 
less, it is more convenient to use a simplified procedure. We can indeed avoid having to 
explicitly integrate over the scattering angle, by first performing the angular integration 
of the leptonic tensor. Doing this, we can compute as and a a by simply performing the 

"Throughout this paper we explicitly consider the case of the 6-quark. The results for the charm quark 
can be simply obtained by properly replacing the mass and the electroweak couplings of the massive quark. 
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integration over the final-state parton momenta of the following projections of the hadronic 
tensor: 



a s oc -g^T^, (6) 
n x Q p 

a A oc ie^xp 77 T pv , (7) 

n ■ Q 

where Q p is the total incoming momentum and the light-like (n 2 = 0) vector n M identifies 
the forward direction. 



2.1 Leading order 

At the leading order (LO) we have to consider the cross sections for the process e + e~ — > bb 
at the tree level and thus, the 6-quark direction and the thrust direction coincide. The 
tree-level cross sections and are straightforward to calculate and the result is[| 



a 



(0) 



rr (0) 



47ra 2 iV, 



,(Q 2 -M 2 )Q 2 



sq 2 



1 e 2 e P v el + 2 y ^ - (P v e e v e + P a e e a e )e b v b 



D Z {Q 2 



Q 4 



D Z (Q 2 

Ana 2 N c 3 \ (Q 2 - M 2 )Q 2 



3Q 2 



4^ 



D Z (Q 2 

Q 4 



v\ + a 2 e )P v + 2P a v e a e (v 2 + a 2 ), 
(P v e e a e + P a e e v e )e b a b 



D Z (Q 2 



2P v v e a e + P a (v 2 + a 2 ) 2v b a b 



(8) 



(9) 



with 



Pv 
Pn 



1 + PlPr , 
Pl + Pr, 



(10) 

(11) 



where Pl is the left-hand-polarization of the electron (+1 = fully left-handed, = unpo- 
larized, -1 = fully right-handed) and Pr is the right-hand-polarization of the positron (+1 
= fully right-handed and so forth), 



D z (Q 2 ) = (Q 2 -M 2 z y + (T z Mz 



(12) 



Ci is the electric charge in units of the proton charge (i.e. e e = —1) and the electroweak 
couplings are: 



1 



2 sin 9 W cos 9 V 
1 

2 sin 6L, cos 6L 



-{t^ - 2e i sm 2 6 w ) , 



t3i 



(13) 
(14) 



T Unlcss explicitly mentioned, we neglect the 6-quark mass throughout this paper. At LO the dominant 
mass corrections are proportional to m 2 / Q 2 and can be found, for instance, in Ref. B . 
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The ratio between Eqs. and (Q) is insensitive to the fine structure constant a and the 
number of colours N c and thus, at LO the forward-backward asymmetry Ap B , 

AO) _^A_ 
^-FB — (o) ' 

gives a direct measurement of the electroweak couplings. In particular, if we are exactly 
on the resonance, Q 2 = Mf , and we neglect the photon contribution, we obtain 

(0) _ 3 Ae + V . , , 

Afb ~ 4 TTA^p Ab ' (16) 

where 

Ai = — 2 ; o > (1^) 

V ~T V - TTp~lPr~ ' (18) 

Finally, for unpolarized beams, we obtain 

4°B = l^eA b . (19) 

This is the form in which the forward-backward asymmetry is most often presented. It is 
worth pointing out however that all of our results will be universal multiplicative correc- 
tionsfj], so apply equally well to any of the forms ( |I5| , [I~6| ) or (|19|). 

Another important variable is the so-called left-right forward-backward asymmetry M, 

A N F (V = +1) - N F (V = -1) - N B (V = +1) + N B (V = -1) 

LR,FB ArF( p = +1) + ArF( p = _ 1)+iVB( p = +1) + iVi?( p = _ 1) • I 

Its LO expression can be obtained from Eqs. (§) and (]£]), and again neglecting the photon 
contribution exactly on the Z resonance, it is given by: 



3 
4 

Our results apply equally well also to this observable 



Af R , FB = \A b . (21) 



2.2 Next-to-leading-order corrections 

At next-to-leading order (NLO), we have to consider the one-loop cross sections (jC 1 ); 011 * 3 "! 00 ? 
for the two-parton process e + e~ — > bb and the tree-level cross sections a^' tree for the three- 
parton process e + e~ — > bbg. We obtain: 

rr^ _i_ ^.(l);one-loop _(l);tree 

4(1) = °A +°A + °A (2?] 

FB 4_ „(l);one-loop (l);tree ' 1 ; 



*At 0(a§) there are some non-universal corrections, but we do not explicitly compute them (see the 



discussion in Sects. pjA and 3.1) 
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Each of the cross sections at 0(«s) is separately divergent, so they have to be regularized 
in some way before being combined together. In any regularization scheme that preserves 
the helicity conservation of massless QCD| (for example, dimensional regularization), we 
have the property 

(l);one-loop (l):one-loop 

= 2 *-?>— • < 23 > 

and hence, if we expand the ratio in Eq. fl2"2"|) up to 0(as), the one-loop corrections cancel, 
and we obtain 

a (0) / o _(l);tree ^.(l^treex 

Although o"^' tree and a s ^' tree are each separately divergent in the soft and collinear regions, 
the divergences cancel at the integrand level, and the whole thing can be calculated in the 
unregularized theory. 

At this order, the different definitions of the forward-backward asymmetry give different 
results. As already anticipated, we consider two possible definitions of the forward direction: 
the 6-quark direction and the thrust axis direction. 



It is straightforward to calculate the NLO corrections in Eq. ( p4[ ) analytically with either 
definition. We obtain: 



4$ = A% (l - \ Cf ^ * Af B (1 - 0.318 a s ) , (25) 

= 4%(l-{j-41n^ + ^ + ln 2 2-^ln3 + 2Li 2 (-i)}^) (26) 

~ AP B (l - 0.670CV^ ~ A% - 0.285 a^j . (27) 

The result in Eq. ( p5|) is well known ||. The analytical result in Eq. (|26"D agrees with the 
numerical calculation performed in Refs. || |6|. The difference between the two definitions 
is only about 0.4% for as ~ 0.12. 

We remind the reader that the NLO QCD correction to the symmetric cross section erg 
in the massless limit is equal to the correction to the e + e~ total cross section, namely 



a s = a%> ( 1 + -C F f + O(o4) ) . (21 



Thus, Eqs. (|25| ) and ( p6| ) imply the following results for the antisymmetric cross sections 
a\ = a 



" ~ J0) (l + 0(4)) , (29) 
a T A = ai O) fl-{l-41nU^ + ln 2 2-^ln3 + 2Li 2 (4))^U0(a|)') (30) 



af (l + 0.034a s + C(ai)) • (31) 



§Note that the relation is explicitly violated for massive quarks. 
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The vanishing of the 0(as)-correction to the antisymmetric cross section a\ with respect 
to the fe-quark axis in the massless case was first noticed in Ref. ||. 

Unlike at LO, the corrections to A^) B due to the finite mass of the b quark are of 
Oim^/Q). The mass corrections have been computed in analytic form for the 6-quark 
direction Q and numerically for the thrust direction ||. 



3 Contributions at next-to-next-to-leading order 



At next-to-next-to-leading order (NNLO) we have to consider the diagrams of Figs. [T]-[5]. 
The single diagram drawn in Fig. |I]b stands for all the one-loop diagrams with one virtual 
gluon. Analogously, the diagram in Fig. [| stands for all the tree-level diagrams contributing 
to the bbgg final state, and so forth. 

We separate the contributions to the cross sections into three classes: flavour non- 
singlet (NS), flavour singlet (S), and interference (or triangle) (TV). We thus write the 
cross sections as 

<?s = °s,ns + 4 2 5 + 4 2 k + 0(<4) , (32) 

OA = VA,NS + 4 2 Tr + 0(<4) ■ (33) 

(2) (2) 

In this notation, up to 0(ats) there are only non-singlet contributions. Thus, &ssi a STr 

and cr^Tr are proportional to «§■ Note also that there are no singlet contributions to the 
antisymmetric cross section a a. 




Figure 1: Some of the diagrams contributing to the bb final state up to C(a|). The dashed 
line represents either the axial or vector current, the thick line the b and the thin line another 
quark q, which must be summed over flavours, including the b- and t-quark contributions. 
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Figure 2: Some of the diagrams contributing to the bbg final state up to 0(a$). The dashed 
line represents either the axial or vector current, the thick line the b and the thin line another 
quark q, which must be summed over flavours, including the b- and t-quark contributions. 




Figure 3: One of the diagrams contributing to the bbgg final state at C(a|). The dashed 
line represents either the axial or vector current. 




Figure 4: Some of the diagrams contributing to the bbqq final state at 0(a^). The dashed 
line represents either the axial or vector current, the thick line the b and the thin line some 
other quark flavour q, with q^b. 




Figure 5: One of the diagrams contributing to the bbbb final state at 0(a$). The dashed 
line represents either the axial or vector current. The cross indicates which of the two b 
quarks is triggered on: both contributions must be summed. 
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The forward-backward asymmetry is decomposed in a similar way. Expanding the ratio 
cta/o's u P to 0(al), we write 

_(0) / (2) (2) (2)\ 

4(2) _ 4 (2) , a A a A,Tr _ <?S£r _ a S,S m \ 

FB - FB,NS -T (o) (0) (0) (0) ' ^ 

a S \ a A a S a S } 

where Ap B NS denotes the non-singlet component: 

A (2) _ (TA,NS , , 

^FB^S — ■ \°°) 

&S,NS 



We now discuss our treatment of each contribution in turn. The classification of the 
four-6 contribution of Fig. |5] also warrants additional discussion. 



3.1 Triangle contributions 



In this class we group all the cross section contributions consisting of two quark triangles, 
one attached to each current. These correspond to the interference between the diagrams 
in Figs. |I|d and [I]a, between those in Figs. ^|c and ^]a, and between those in Figs. [|b and [|a. 
They give non-universal (i.e. non-factorizable) corrections to both the symmetric and anti- 
symmetric cross sections. They are calculated in Ref. |7[] for the 6-quark axis definition and 
found to be very small. To our knowledge their contribution to the thrust axis definition 
has never been calculated, but we expect it to be similarly small. We therefore neglect it, 
i.e. <Jg Tr and cr^ Tr in Eq. from our calculation^]. 



3.2 Singlet contributions 



In this class we group the square of the diagrams of Fig. ^b, where the final-state b quark is 
not coupled to the current. In these contributions the b and b are produced in a definite state 
of charge conjugation, C = +1. They therefore cannot contribute to the antisymmetric 
cross section, cta- Their contribution to the symmetric cross section, as, is logarithmically 
enhanced in the small-mass limit and proportional to agin 3 Q 2 /m^. An approximate ex- 
pression for it, denoted by F BTanco , was used in Ref. 0. It is calculated exactly to 0(a\) 
in Refs. ||T0| , pf , and the leading and next-to- leading logarithms are summed to all orders 
in as in Ref. [ |TT|j . 

Note that the singlet contributions to a,s include an additional term coming from the 
bbbb final state. As discussed in Sect. |373|, this term is very similar to that described above. 
It was missing in the expression denoted by F Branco in Ref. 0. 

In some sense the singlet component is a 'background' to the forward-backward asym- 
metry measurement and, in fact, in the experimental analyses (see e.g. Ref. |T2|]) it is 
statistically subtracted using Monte Carlo event generators. We therefore neglect it, i.e. 
a ss in Eq. ([34]), from our calculation. 

^We remind the reader that the triangle contributions to both us and <ta are finite in the massless 
limit mj — > 0, provided that the sums over quark flavour q in the diagrams of Figs, ^d and run over 
complete SU(2) doublets. 
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3.3 Four-6 contributions 



The classification of the four-5 diagrams of Fig. |5] deserves special mention. Let us first 
point out a basic fact. The four-6 diagrams of Fig. |5| contribute to both the 6-quark 
cross sections as and a a and the e + e~ total cross section. However, they appear with 
different multiplicity factors in the two cases. In the case of the e + e~ total cross section 
the multiplicity factor is simply equal to unity. In the contribution to the inclusive 6-quark 
cross sections 05 and a a, these diagrams count twice since there are two b quarks in the 
final state. This observation is relevant in the discussion that follows and, in particular, it 
is important in understanding the results for the non-singlet component of the symmetric 



cross section as discussed in Sect. 3.4 



After summing and squaring the diagrams in Fig. [5], we obtain two types of contribution: 
i) those that are identical to the contributions of Fig. [| but with the other quark q replaced 
by an untriggered-on b quark, and ii) those that are genuine interference terms arising from 
the fact that the two antiquarks are indistinguishable, called the E-teim. in Ref. |13|]. The 
squared diagrams of type i) are treated as those of Fig. f|, that is, we lump them together 
with the corresponding terms from Fig. |] in the sing let (afs in Eq. ©), non-singlet 



{&s,ns and aA,NS in Eqs. (|32|) and (|33|) ) or triangle (crlf^y and a^ Tr in Eqs. (|32"1) and (|33|) ) 
contributions. The squared diagrams of type ii), which give a universal (i.e. factorizable) 
correction to both the antisymmetric and symmetric cross sections, can be considered part 
of the non-singlet contributions. 

It is not entirely clear how four-quark final states are actually treated in the different 
experimental analyses, i.e. the extent to which they are genuinely measuring the inclusive 
cross sections. Often some vague statement like "a four-6 final state is more likely to 
be tagged than a two-6 one, but less than twice as likely" is made. To know what to 
calculate one must understand the corrections that are applied for this difference in tagging 
efficiency, which are not usually explicitly stated in the papers. In the absence of a unique 
experimental procedure and of a definitive statement from the experiments on what they 
are measuring, we make this ambiguity explicit by multiplying the E'-term by an arbitrary 
weight factor We 0- An inclusive definition would correspond to We = 2 (each b quark 
contributing once), while an exclusive definition (the cross section for events containing at 
least one b quark) would correspond to We = 1- Since the forward-backward asymmetry 
is defined to be the asymmetry of a differential cross section, it is clear that we must use 
the same cross section definition in the numerator and denominator, i.e. that We must be 
the same in the symmetric and antisymmetric cross sections. 

We return to the role of the weight factor We after discussing the general form of the 
non-singlet contributions. 



"Note that we use the same normalization as in Ref. ]13| (see also Eq. (BJ)) in which the E-terni 
already includes an identical-particle factor of 1/(2! ) 2 because there are two identical quarks and two 
identical antiquarks in the final state. Thus, when we set We = 2 we actually include an overall factor of 
W E /(2\f = 1/2!. 
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3.4 Non-singlet contributions 



Here we consider all the other contributions that have not yet been treated, namely all the 
diagrams in Fig. [I] except those in Fig. |l]d, the diagrams in Figs. 0a, 0b, |3|, [|a and ||a, as well 
as the E-term defined above. All these terms are included in the non-singlet components 
&s,ns and &a,ns of Eqs. (^) and (|33]). Actually, introducing the weight factor W F for the 
E-term, we can define the following symmetric and antisymmetric cross sections 

<y S ,Ns{W E ) = <ts,ns{We = 0) + W E 4° } / E s , (36) 

va,ns{We) = a ANS (W E = 0) + W E af J E A , (37) 

where / Es and / E A denote the integral of the symmetric and antisymmetric E-term, 
respectively. We recall that the 'truly' inclusive cross sections in Eq. (f|) correspond to the 
definition with We = 2, i.e. <ts,ns — &s,ns(We = 2) and (Ja,ns = <?a,ns(We = 2). 

The (9(ag)-calculation of the cross sections in Eqs. (|36| , [37]) and of the corresponding 
forward-backward asymmetry in the case of a finite 6-quark mass is extremely complicated, 
and we are not able to perform it. It is thus convenient to separate the calculation into a 
piece that is finite in the massless limit and a simpler piece that is not. Then, the (although, 
cumbersome) finite piece can be more easily computed in the massless approximation, while 
the simpler non-finite piece can be computed in the massive theory. 

It is possible to show (Appendix A) that the inclusive definition, with We = 2, results 
in an antisymmetric cross section oa (or, analogously, <ta,ns) that is finite in the massless 
limit, at least at 0(a|). However, in the same limit, the inclusive symmetric cross section 
is divergent at (9(a§), even if we only consider its non-singlet component. The corrections 
to (the non-singlet component of) the forward-backward asymmetry itself must therefore 
also be divergent in the massless limit. 

This final statement remains true for any value of We > 0. For example, with We = 1, 
the non-singlet part of the symmetric cross section is finite (see Eq. ([¥!])), but the antisym- 
metric cross section contains logarithmically-enhanced terms. 

The divergences in the non-singlet components correspond to logarithmically-enhanced 
terms a\\nQ 2 jm\ coming from the .E-term in the triple-collinear limit, i.e. when three 
fermions of the four-quark final state become simultaneously parallel (Appendix B). The 
integral of the symmetric .E-term is calculated for finite values of the quark mass in Ap- 
pendix B. Neglecting corrections of 0(rrib/Q), the final result is 



/ 



2 




2 



2C 3 j In % - 8.1790 ±0.0013 



rn 



(38) 



As expected from the singular behaviour in the triple-collinear limit, the analytic coeffi- 
cient in front of lnQ 2 /ml is proportional to the integral of the non-singlet Altarelli-Parisi 
probability Pjl s (z, as) (see, for instance, the first paper in Ref. 



m 

\zP^{z,as) = (^) 2 C F (c F - 1 -C A ){ 1 4-^ + 2C 3 ) (39) 



o 



2 




^) C F [C F -^-C A ) 0.7193 . (10) 
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The constant term in the square bracket on the right-hand side of Eq. (|38D is the result of 
our numerical calculation. 

Having pointed out that the symmetric i?-term is divergent in the massless limit, it is 
very simple to show how the divergence appears in the inclusive symmetric cross section. 
According to the definition of the non-singlet component of as, the virtual diagrams that 
contribute to cts,ns are exactly those that contribute to the e + e~ total cross section. As 
for the real diagrams, they only differ by the contributions coming from the i?-term. In the 
total cross section, the i?-term enters with a multiplicity factor We = 1, and its divergence 
is cancelled by that of the virtual diagrams. In the inclusive 6-quark cross section, the 
multiplicity factor of the .E-term is We = 2 and, thus, the cancellation of the divergence 
with the virtual terms is spoiled. 



This argument also allows us to directly compute the C(ag)-correction to Eq. 
Exploiting the fact that the massless QCD correction to cts,ns(We = 1) is equal to the 
correction R e + e - to the total cross section, we write 



&s,ns(We = I) = (t ( s 



(0) 



R P 



(41) 



and, more generally, 



<ts,ns(We) = o 



(o) 



Re+e- + {We - 1) / E s + 0{a% 



(42) 



Then, we obtain an explicit expression for <ts,ns(We) by simply introducing in Eq. 
our result in Eq. fl38l) for J Es and the well-known result [15| for R e + e -'- 



R» 



3 a s (Q 2 ) 
1 + 4° F —^ 



+ 



' as(Q 2 
2tt 



iv r 3 

-j C f I—-C f + Ca 



where Tr = 1/2 and Nf is the number of light flavours (e.g. Nf = 5 at LEP) 
In particular, for the inclusive symmetric cross section we obtain 



&S,NS = &s,ns(W e = 2) = a s 



(0) 



Re+e- +J E S + 0(<4) 



(43) 



(44) 



The explicit C(ag)-calculation of the antisymmetric cross section <ja,ns and of the 
forward-backward asymmetry is described in the next Section. 

Note that our result in Eq. ((44]) for the inclusive symmetric cross section disagrees 
with the corresponding result of Ravindran and van Neerven M. Their expression for the 
correction to the symmetric cross section (fy + fi in their Eqs. (31) and (32)) is equal to 
the result in Eq. ( |4"3"| ) for the C(a|)-correction to R e + e -. The disagreement thus regards the 
additional logarithmically-enhanced term / Es included in our expression. The multiplicity 
of 6-quarks is not required to be finite in massless QCD (even in the non-singlet sector), 
and thus we cannot find any reason why this logarithmically-enhanced term can disappear 
from the inclusive symmetric cross section. 
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The results of Ref. || for <ts,ns — &s,ns(We = 2) are based on the calculation of the 
single-particle inclusive distribution performed in Refs. fll6| . Our result is consistent with 
those in Refs. ||16|| . In fact, we have evaluated the integral over the longitudinal-momentum 



fraction z of the non-singlet coefficient function Cg^(z,as(Q 2 ),Q 2 /^ 2 F ) = C^q + Cf & 



y L,q i 

computed there. This integral is proportional to (Js,ns — &s,ns(We = 2) in massless QCD 
after factorization of collinearly-divergent contributions at the factorization scale [ip. We 
find that the integral explicitly depends on lnQ 2 /n F , thus proving that (Js,ns is n °t finite 
in massless QCD. The coefficient of In Q 2 / fi F exactly agrees with the coefficient of In Q 2 /m 2 
in our Eqs. (p8T) and (E3). 



4 Calculation of the non-singlet contribution at 0(a|) 



As discussed in Sect. |3.4| , the NNLO corrections to the non-singlet component of the 
forward-backward asymmetry, Afb,ns, are divergent in the massless limit. The diver- 
gent behaviour remains true also if we abandon the fully inclusive definition by introducing 
the arbitrary weight We for the E'-term. Thus, A F b,ns cannot be computed at 0(a\) by 
using the massless approximation. 

Nonetheless, since both cta,ns(We — 2) and (Js,ns(We = 1) are finite when m b — > 0, 
we can use the dependence on We to construct an unphysical observable that is finite in 
the massless limit: 

^(2);finite _ ct a ,ns(We = 2) 
FB cr s ,Ns{W E = 1) 

The physical result for We = 2 is then given by 

4 (2) _ ,(2);finite a a* f p _ 4 (2);finite .(0) f p (Ap s 
A FB,NS - A FB ~ J hs - FB ~ A FB J , (4b) 



S 



where / Eg is the integral of the symmetric E-tevm, given in Eq. (^). 

The massless calculation of A FB ^ n ^ e can be performed in a similar way to the NLO 



calculation of Sect. 2.2. The total contribution can be written as 



. _(0) (1) (2) ; two-loop (2) ; one-loop (2) ; tree m/ _ 9 \ 
.(2);finite _ O A + a A + a A + a A + ° A [W E - Z) 

FB " _(0) , _(l) , (2) ; two-loop (2) ; one-loop (2);tree/ u/ _ ,\ ' { ) 

a S + a S + a S + a S + a S \ W E — J- J 

where and are the complete contributions to the antisymmetric and symmetric 
cross sections at 0{ots). The non-singlet (9(ag)-contributions from the two-parton, three- 
parton and four-parton final states are denoted by (j(. 2 ); two -\ 00 P ; a (2) ; one-loop anc j CT (2);tree 
respectively. Of course, the dependence on We enters only through the four-parton terms 
af' txee (W E = 2) and a^ tiee (W E = 1). 

If we continue to use a regularization scheme that preserves the helicity conservation of 
massless QCD, like dimensional regularization, the two-loop corrections are again propor- 
tional to the tree-level results, 

(2);tW0-l00p (2);tw0-l00p 

ZA - (A9\ 

JO) (0) ' 

°A a S 
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so that if we expand the ratio in Eq. (f47|) up to (9(a§), the two-loop corrections cancel, 
and we obtain 



.(2) ; finite _ a A 
A FB — — 



(0) 



a 



(0) 




a 



s 



(2);one-loop 
a A 



(2);one-loop 
a S 



a 



.(0) 



a 



.(o) 



+ 



of) 

af' tree (W E 



2) 



a 



(2);tree 
S 



a 



(0) 



a 



(0) 



.(49) 



The first line can be calculated analytically (see Sect. [2.2| ), but the second line is too 
complicated to be able to, so must be done numerically. Since the two-loop terms have 
cancelled, this has the structure of a NLO three-jet calculation, as first noticed by Altarelli 
and Lampe JF]. Thus the calculation can be performed using known techniques (we use the 
dipole-formalism version of the subtraction method |]17|| ). One simply has to replace the 
full matrix element squared by the appropriate contractions of the hadronic tensor, as in 
Eqs. (|6|) and (0). We have obtained simplified analytical expressions for these contractions 
by using the matrix elements originally computed by the Leiden group 
checked that these expressions numerically agree with the code of Ref. 



T8|. We have also 
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4.1 Numerical results 



We are finally ready to present our numerical results. We start with the unphysical, but 
finite, quantity defined in Eq. fl45|) , and separate out the different colour factors, as in 
Refs. 



. (2);fmite;6 
A FB 



A 



(0) 
FB 



1 _ «S 
2vr 



«s 3 

u p 

2tt2 



3 -c F 

2 



+ (g) 2 C F (CC F + NN C + TT R Nf 



(50) 



with as = asiQ 2 )- Our numerical results are shown in Table [I], in comparison with the 
previous calculations. It is clear that we disagree badly with the results of Altarelli and 
Lampe M, but are in excellent agreement with Ravindran and van Neerven M, who give 
the coefficients analytically. However, we should recall that this must have subtracted from 
it the logarithmically-enhanced term of Eqs. fl4"E| , 
of Ref. 



which is not present in the result 



In fact, in Sect. |3.4| we have already pointed out that their expression for the 
correction to the symmetric cross section does not agree with ours, but, rather, it is actually 



6-quark axis 


C 


N 


T 


AL i 


4.4 ±0.5 


-10.3 ±0.3 


5.68 ±0.04 


RvN § 


| = 0.375 


-If. = -15.375 


T = 5-5 


Our Calculation 


0.3765 ±0.0038 


-15.3769 ±0.0034 


5.5002 ±0.0008 



Table 1: Results for the coefficients of the O(o§) correction to the finite part of the forward- 
backward asymmetry with the 6-quark axis definition, Eqs. (|4*6| , |50| ) . 
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thrust axis 


C 


N 


T 


Our Calculation 


-3.7212 ±0.0065 


-9.6011 ±0.0049 


4.4144 ±0.0006 



Table 2: Results for the coefficients of the 0(a\) correction to the finite part of the forward- 
backward asymmetry with the thrust axis definition, Eqs. (|46| , |5l|) . 



equal to our 



(We = 1)- So, the fact that their result for the complete Ap B agrees 
with our v4^^ nite means that we confirm their result |20], |8| for the inclusive antisymmetric 



cross section = <t%'[We = 2) (f A in Eq. (33) of Ref. §). 

The disagreement with the result of Ref. J/J may be related to the poor numerical 
convergence of their calculational method (i.e. the effect of large numerical cancellations). 

Using our numerical program it is straightforward to calculate the forward-backward 
asymmetry with any other axis definition (or cuts, for example on the value of the thrust). 
With the thrust axis definition, we obtain 



r( 2 ) 



.(2);finite;T _ AO) 
s±FB — ^FB 



as V 



( ^ J C F (CC F + NN C + TT R N f ) 



(51) 



with as = ctsiQ 2 ) and the coefficients given in Table |2|. The logarithmically-enhanced 
piece that has to be added to this is identical to that in the 6-quark axis definition, namely 
Eqs. ( [46"| , |38|) . It is worth noting that the difference between the two definitions is the same 
size and in the same direction as at 0(a>s), leading to an overall difference of 0.8% for 
a s ~ 0.12. 

Since Ap^ n ^ te ' T is defined by the ratio in Eq. (|45|), using the expression in Eq. ( f4~T|) for 
vs,ns(We = 1), we can translate our result in Eq. ([y]) into an equivalent result for the 
antisymmetric cross section defined with respect to the thrust axis. We have: 



a 



A,NS 




11 

2 



+ (^-ll( 3 + N)C A +U( 3 -^+T)T R N f 



with the coefficients C, N and T given in Table 0. 




(52) 



We finally recall that we include an arbitrary factor We in front of the four-6 contri- 
bution to account for the way in which it is treated in the experimental analyses. For a 
fully inclusive definition, in which each b quark contributes once, We should be set equal 
to 2, while for an exclusive definition, We should be set equal to 1. Our final result for the 



**In the analogous expression for <J ANS , the coefficient of as(Q 2 ) vanishes and C,N and T are those 
given in Table [l], which exactly cancel the rational numbers in Eq. (|5^), leaving only 3/3qC3j with /3q = 



TrNj, as pointed out in Ref. J2C 



14 



non-singlet component of the forward-backward asymmetry, is then: 



4 (2) (w x _ <7A,Ns(We) _ 4 (2);fimte ,(0) 
A FB,NS\ VV E) = — /yi7 \ — A FB ~~ A FB 



(1 



\W E ) 



2JEa- J e s)+ \W e J E s 



(53) 

where A^p B mte is given in Eqs. ([5(], [51]) and Tables [1] and [| / £s is given in Eq. ( |3~8"D and 
(see Appendix B) 



2jE A - Je s = (g) C F (C F - iC A ) (0.3620 ± 0. 



0007 



2 E A ~ E 



^Yc F (C F - ±CV) (0.1144 ±0.0009 



Note that the combinations of .E-term contributions in Eqs. (|54[) and 
massless limit (see the discussion in Appendix B). 



quark axis, (54) 
thrust axis. (55) 
5|) are finite in the 



Putting all these numbers together, and setting Nf = 5, we write the forward-backward 
asymmetry according to the two definitions as: 



A (o) 

FB 



A 



(2);T 
FB,NS 



{W E 



A 



(o) 

FB 



1 - 0.318a s - 0.973c4 + W E a 2 s [0.00405 In ^ ~ °- 0240 

V m b J 

1 - 0.284a s - 0.676a^ + W E a\ (0.00405 In % - 0.0233 



, (56) 
• (57) 



5 Conclusion 



We have calculated the second-order corrections to the non-singlet component of the 
forward-backward asymmetry in e + e~ annihilation. We have retained all terms that do 
not vanish in the small-mass limit (constants and logarithmically-enhanced terms). Our 
result is also valid for the left-right forward-backward asymmetry. 

Using the quark axis definition we do not agree with any existing calculation. Separating 
the asymmetry into its symmetric and antisymmetric parts, we find that we agree with 
Ravindran and van Neerven || for the antisymmetric part, which is finite in massless 
QCD. For the symmetric part we disagree by a term that is divergent in massless QCD, 
so is logarithmically-enhanced in the full theory. 

We have obtained results for the first time with the thrust axis definition, which is 
actually what is used in the experimental measurements. After including the second-order 
contributions, the difference between the two axis definitions is twice as large as at first 
order, amounting to 0.8%. 

We summarize the total QCD correction according to the various available calculations 
in Table [3[ We continue to neglect all terms that vanish in the massless limit, and discuss 
the effect of mass corrections below. Since in the existing experimental analyses (see for 
example Ref. ||12|)) the known 0(ats) correction for the thrust axis definition was included, 
together with the Altarelli and Lampe quark axis value for the O(ot\) corrections, we do 
the same in Table |3|. 
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We find that the difference between the Ravindran and van Neerven calculation and ours 
is numerically irrelevant, being smaller than 10 -4 for b quarks and ~ 2.5xl0 -4 for c quarks. 
Therefore at the numerical precision required by current or any foreseen experiments, we 
agree with their result - the difference is only one of principle. The difference between the 
Altarelli and Lampe calculation and ours for the quark axis definition is more significant 
though, at around 1%. However, the error in their calculation and the effect of using the 
thrust axis definition partially cancel, and the total difference is around 0.6%. 

Before quantifying the impact of these differences, we mention the important fact, 
discussed in Ref. [|12|], that the experimental procedures introduce a bias towards more 
two-jet-like events. This actually decreases the size of the QCD corrections considerably, 
so our numbers should be considered as upper bounds. In fact at present the effect of 
this bias is typically taken into account using Monte Carlo event generators. Using our 
numerical calculation, it is straightforward to apply any infrared-safe cuts, for example on 
the thrust of the event (this effect was first considered at 0(as) in Ref. ||). This could be 
used to reduce the reliance on the generators, or at the least to calibrate their reliability. 

To quantify the impact of the differences shown in Table |3|, we recall a few figures from 
the latest global electroweak fit 0. The total error on the LEP average forward-backward 
asymmetry of 6-quarks Ap B is 2.1%. The second-order QCD corrections are used to convert 
the measured value into a measurement of the tree- level asymmetry, Ap B , at present using 
the Altarelli and Lampe result. This is then used as input into the fit for the effective weak 
mixing angle, sin 2 9 e ff, and eventually into the global fit to all electroweak data. Following 
through this process, our smaller value of the correction in Table [| results in a larger value 
of Ap B and hence a smaller value of sin 2 9 e jf, by about a third of its experimental error. 

This has a direct bearing on the fitted value of the Higgs mass, (see Fig. 9 of Ref. fLj). 
To find the effect of using our corrections would require a complete refitting of all the 
electroweak data. However, we can get a rough idea simply by fitting the data in Fig. 9 of 
Ref. |]J alone. We find a roughly linear relation: for each per mille that the corrected value 
of the quark asymmetries is increased, we obtain a per cent decrease in the central value 
of the Higgs mass (and its upper bound). Therefore with our 0.6% difference, we expect 
a reduction of about 5 GeV in the central value. While this is certainly not statistically 
significant, given the importance that some people attach to this value, it is not irrelevant 
either. 

In trying to estimate the remaining uncertainties in the forward-backward asymmetry, 
we recall the ingredients still missing from our analysis. We should bear in mind that while 
the 2% precision of current experiments is close to their final limit, a future linear collider 





AL § 

quark axis 


RvN § 
quark axis 


Our Calculation 
quark axis 


Our Calculation 
thrust axis 


Correction, Ap B /Ap B 


0.962 


0.952 


0.952 


0.956 



Table 3: Total QCD correction to the forward-backward asymmetry in the small-mass limit, 
with as = 0.12. In each case, the thrust axis definition is used for the 0(as) correction 
and the definition shown is used for the 0(a^) correction, as discussed in the text. 
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could be capable of experimental errors on the left-right forward-backward asymmetry of 
order 0.1% §. 

Within small-mass perturbation theory, the first terms that we neglect are 0(a$). To 
estimate their size, we assume that the coefficient grows as much in going from 0(atg) 
to 0(ag) as it did from 0(as) to (9(a|), and get 0.3%. The more conventional method, 
varying the renormalization scale from Q/2 to 2Q results in a similar estimate of 0.2%. A 
variation in the input value of as of ±0.004 gives only 0.1%. 

Within the 0{a\) calculation, we neglected the effect of triangle diagrams. For the 
quark axis definition, these were calculated in Ref. @, and amount to about 0.1%. We 
have no reason to suppose they would be larger for the thrust axis definition, and in any 
case it would not be difficult to calculate them. 

We have also neglected linear mass corrections of the type rrib/Q, which are absent at 
tree level, but arise at higher orders. The full mass correction at 0(as) is well known, and 
is reasonably well approximated by its leading term, 4Cpas/Tim b /Q. Since we do not have 
any higher order corrections to this linear mass term, its renormalization group dependence 
is not under control, so to estimate the effect of higher order corrections, we vary m& from 
its running value in the MS scheme (~ 3 GeV) to its pole value (~ 5 GeV), resulting in a 
0.4% variation in Aps- 

Finally, at higher orders it is quite possible that the leading mass term could become log- 
arithmically enhanced, for instance, as ~ a^rrib/Q ln n (Q 2 /mf) at the second order. Terms 
like this certainly arise with n=l simply from the renormalization group effects just men- 
tioned, but the question is whether additional terms can arise from other dynamic effects. 
A possible additional source of single-logarithmic enhancement is collinear emission, as in 
the case of the E-term contributions discussed earlier. Owing to the inclusiveness of the 
forward-backward asymmetry with respect to soft emission, we think that higher powers 
of logs are unlikely to be present in the non-singlet component at 0(a^). Although this 
point deserves further investigation, assuming n < 1 we estimate a resulting uncertainty of 
0.5%. 

We have not made any attempt to estimate the uncertainty due to non-perturbative 
corrections. In Ref. [TjJ, this is done using Monte Carlo event generators. They find a 
correction of 0.25% and conservatively assign the whole of this as a systematic error. 

To summarize, there are several sources of uncertainty that all contribute at the few per 
mille level. While this is certainly sufficient for the current precision of the data, matching 
the precision of a future linear collider measurement could be extremely difficult. It is likely 
that this could only be done by making even more stringent two-jet cuts in order to work 
in a region in which the corrections and their uncertainties are smaller. 
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Appendix A: The antisymmetric cross section in massless QCD 



In this Appendix we show that, up to 0(a|), the perturbative-QCD corrections to the 
heavy-quark antisymmetric cross section a a are finite in the limit of vanishing quark masses. 

We are interested in the analogue of the cross section in Eq. (|J) for the inclusive process 
e + e~ — ► a + X where a = qf, qf, g denotes a generic massless QCD parton. We thus define 
the antisymmetric^ cross section da A /dx as follows 

da A f 1 j n da ( e+e ~ a + X ) f° j n d<r(e+e- -> a + X) . , 



r^U = x"' 1 ^ , (A.2) 



dx 7o dx dcos 9 J-i dx dcos 9 

It is also convenient to introduce the iV-moments a A N defined by 

- 1 -..dal 
dx 

and likewise for any other function of the energy fraction x. Note that the massless limit 
of the 6-quark antisymmetric cross section in Eq. (|j) coincides with the N = 1 moment of 
da q //dx, i.e. o A = a% N=1 . 

Since we are working in massless QCD, the antisymmetric cross section da A /dx is not 
finite in perturbation theory and, more precisely, it is collinear divergent. Nonetheless, 
because of the factorization theorem of mass singularities, once the divergences have been 
regularized (by using, for instance, dimensional regularization) they can be factorized. The 
TV-moments can be written as 

a A,N = °A,N ^ba,N , (A. 3) 

b=q f ,q f ,g 

where a\ N is a finite contribution to the cross section and the factor T a b t N contains all the 
collinear singularities (see e.g. Ref. [pj}|). This factor depends on the factorization (or reg- 
ularization) scale \i and the factorization scheme, but it is universal (process independent). 
Moreover, it fulfils the Altarelli-Parisi evolution equations 

d 

^ab,N = 5I jP «c,v(as(/U 2 )) Fcb,N , (A.4) 



with the initial condition r^iv^ 2 = 0) = 5 a b and where -P ac ,7v(«s) are the iV-moments of 
the Altarelli-Parisi probabilities, whose power series expansion in as can be computed at 
any perturbative order. 

Note that the antisymmetric cross section a A N is an odd quantity under charge con- 
jugation. Thus we have a A N = —<J A N and, in particular, a AN = —u AN and a AN = 0. 
Analogous relations are valid for a AN . 

We can now consider in detail the massless limit of the 6-quark antisymmetric cross 
section a A , that is, the first moment a AN=1 . Since a AN is C-odd, Eq. ( |A.3| ) gives 



<7%N=1 - ^2^A,N=1 r /'/,7V=l , (A. 5) 



^Exactly analogous definitions hold for the symmetric cross section da^/dx and its N- moments ag 



S.N- 



18 



where 

r /'/,jv=i = r 9/ ,g /)JV =i - Tg fiqf)N=1 . (A. 6) 

Using Eq. ( |A.4|) and the property P ac {&s) — Pac( a s)i which follows from the charge- 



conjugation invariance of QCD, we obtain the following evolution equation for the singular 
collinear factor on the right-hand side of Eq. (|A.5|) 







d In fi- 



i>/,jv=i = [ p q f ,q f „,N=i(as(v< 2 )) ~ Pq fiqfll ,N=i{as{n 2 ))] V f» f, N=1 , (A.7) 



Note that the combination of first moments of the Altarelli-Parisi probabilities on the 
right-hand side of Eq. ( |A.7|) vanishes up to (9(ctg): 



P qfiqf „,N=i(a s (fi 2 )) - Pg fiqfll ,N=i(a s (fx 2 )) = 0{a 3 s ) . (A.t 



This result follows from fermion-number conservation and it can be explicitly checked by 
using the known LO and NLO expressions ]nj of the Altarelli-Parisi probabilities. Equa- 
tion ( |A.8| ) implies that Tpp^i = dp /+0(o!g) and, thus, the massless-quark antisymmetric 
cross section cr q / N=1 is free from collinear singularities up to NNLO accuracy: 

°%n=i = + 0(a 3 s ) . (A.9) 

To conclude our argument on the finiteness of cr^ / JV=1 , we have to discuss the effect of 
soft singularities. The QCD factorization theorem guarantees that the short- distance cross 
section cr q J ( N is finite for any value of the moment index N > 1. The limit N — > 1 of 
Eqs. ( |A.2| ) and (|A.3| ) has to be dealt with with care because it is sensitive to possible soft 



singularities of the type &a q j[ /dx ~ l/x in the inclusive quark spectrum. These singularities 
can arise when qj is produced by the fragmentation of a soft gluon. At 0(as) there are no 
such fragmentation subprocesses. At there is only the subprocess g — qf(jf. In this 

0(ag)-subprocess, however, the pair qfcjf is produced in a definite state of positive charge 
conjugation and, thus, it gives a vanishing contribution to the C-odd cross section a 1 /. It 
follows that up to O(ttg) the limit — > 1 can be safely performed and the right-hand side 
of Eq. flOD is finite. 

The antisymmetric cross section a a is the integral of a single-particle (parton) inclu- 
sive distribution and thus, the finiteness of Eq. ( |A.9| ) may appear surprising. However, 



this result is not accidental. The collinear safety of a 9 / N=1 follows from fermion-number 
conservation, i.e. Eq. ( |A.8| ), and this is exactly the same equation that, up to ©(ctg), guar- 



antees the finiteness of the Adler, Gross-Llewellyn-Smith and unpolarized-Bjorken sums in 
Deep-Inelastic-Scattering. 



Appendix B: Integrating the E-terms 



In this Appendix we give some details of the integration of the i?-terms appearing in, for 
example, Eq. (p3|). 

We begin with the symmetric term / E$, which is equal to the total contribution to 
R e + e - from four-6 final states. We label the quark (antiquark) momenta as pi^ (^3,4) arid 
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retain all mass terms. The integral is then analogous to Eq. (B.2) of Ref. [ 13 1 , but with 
the massless phase space replaced by that with four equal final-state masses: 



(2!) 2 ~ r vS) <^2 I rf$ 4(Q 2 ;^,m^,m^m^ 



Es(pi,P2,P3,Pi) + (pi <-> P2) + (pa <-> Pa) + (pi <-> P2,P3 *-> Pa) 



(B.1) 



and with the E-term itself in Eq. (B.7) of Ref. [HJ replaced by: 



Ec 



(C F -C A /2)< 



Sl2S23 s 34 ~ S12S24S34 + S12S14S34 + S12S13S34 + Si3S 24 



-S14S23S24 + S13S23S24 + S13S14S24 + S13S24 - S14S23 _ S 14 S 23 - S13S14S23) 

— 2m^(2si2Si3 + 3S12S14 + S12S23 — S12S24 — S12S34 + 2si 3 Si4 + S13S24 

+ S 2 13 + S14S23 + S14S24 + S14S34 - S23S24 - S23S34 - S23 - 3^24^34 - S34) 

+4ml(s 12 + 2s 13 + 5s M + s 23 - 2s 24 - 3s 34 ) 



-16m b J/(si 3 S23Sl23Sl34) ~ ^i2(si 2 S 3 4 - S23S24 ~ S13S24 ~ S14S23 ~ S14S13) 

-2m^(4si2S34 + 2s\ 2 - 2s 13 s 14 - 2s 13 s 23 - 2s 13 s 2 4 - s? 3 ~ 2si 4 s 2 3 - 2s 23 s 2 4 - s^) 
-Aml{-Qs 12 + 2s 13 + s u + 2s 23 + «24 - 3s 34 ) 

-24mf V(Si3S 2 3S? 23 ) - ({S U + S1 3 )(S24 + S23)S34 

-ml(2s 12 s 13 + 2si 2 s 23 + S12S34 + s\ 2 + 4si 3 s 2 3 + Si 3 s 2 4 

+ S13S34 + + S14S23 — S14S34 + S23S34 + S23 — S24S34) 

-2m A b (-2s l2 - si 3 + 3su - s 2 3 + 3s 2 4 + 2s 34 ) 
+24m^ J / (S13S23S134S234) 



+ 



(pi ^P3,P2 ^Pa) 



(B.2) 



where Sy = (pi + Pj) 2 and Sijk = (p% + Pj + Pk) 2 — rn 2 . Note that by setting mj = we 
trivially recover the result of Ref. ||13||. 



In the massless case, the integral is divergent in all four triple-collinear limits. When i, j 
and k are all collinear, we have s^- ~ Sik ~ Sjk ~ s^fe — > and the leading behaviour of the 
squared matrix element is ~ \ j s 2 ^ k . Since the volume of three-body phase space is ~ Syk, 
we obtain a logarithmic divergence. Its coefficient is the integral of the corresponding 
Altarelli-Parisi splitting function (either P^ s or P^ q s , which are equal because of the 
charge-conjugation invariance of QCD). After summing over the four singular regions, we 
obtain one singular contribution for each of the two partons in the tree-level contribution, 
so we expect the coefficient of the logarithmically-enhanced term in Eq. ( |B.1| ) to be 



dzP, 



NS( 
qq \ 



(B.3) 
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Figure 6: The term c in Eq. ( |B.4| ) as a function of mass. The errors are purely from Monte 
Carlo statistics. The solid line is our fit to the limiting value and the dotted lines its error. 



That is, we expect the result retaining the quark mass to be of the form 



J E s = C F (C F -C A /2) 



as 
2tt 



13 7T 2 \ Q 2 

2 — -— + 2C 3 ln^ + c 



(B.4) 



with c tending to a constant at small masses. Our numerical results confirm the coefficient 
of the log. For the constant term we obtain the results shown in Fig. |6|. To obtain the 
limiting value, we have tried fitting various degree polynomials in m^/Q to the points 
iTib/Q < /i max , reducing /i max until the fit is acceptable. We call the range of values from 
the different fits a systematic error, which is comparable to the statistical error, and add 
them in quadrature, to give: 

c= -8.7190 ±0.0013. (B.5) 

At the Z peak, logs of the bottom quark mass are not yet asymptotic. Using the 
log-plus-constant approximation, we obtain 

J E s (m b /Q = 5/91) w C F (C F - C A /2) (^j * [-0.3719 ± 0.0015] , (B.6) 

while direct integration gives 

J E s {m b /Q = 5/91) = C F {C F - C A /2) (^j ' [+0.8174 ± 0.0001] . (B.7) 

Even so, the difference between the two results is still an order of magnitude smaller than 
a\m b /M Z) the anticipated size of mass corrections. 
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For charm quarks however, the log-plus-constant approximation works quite well: 



J E s {m c /Q = 1.5/91) « C F {C F - C A /2) [3.0922 ± 0.0015] , (B.8) 

J E s {m c /Q = 1.5/91) = C F {C F - C A /2) Q^V [3.3527 ± 0.0003] . (B.9) 

We turn now to the integral (2 / E A — J Es), which we claim is finite in massless QCD. 
If it is defined in the most natural way, Eq. ( B.10| ), the integrand is not piece- wise finite, 
making it unsuitable for numerical integration. However, we can rewrite it in a form in 
which it is, proving the finiteness of the whole integral, and allowing it to be performed 
numerically. 

If we define E A (n) to be the i?-term contribution that is antisymmetric with respect to 
the direction n, then our integral for the quark axis definition is 

J (E A ( Pl ) + E A {p 2 ) - E s ) . (B.10) 

In each of the four triple-collinear limits — > 0, the integrand diverges like l/sfj k , again 
yielding a logarithmic divergence. The coefficient of this divergence is either positive or 
negative, depending on whether the collinear partons ijk are qqq or qqq. 

However, using the fact that E A is C-odd, we have the relation 

J E A ( Pl ) = J E A (p 2 ) = - J E A (p 3 ) = - J E A (p 4 ), (B.ll) 



which we can exploit to rewrite Eq. ( B.10|) as 



J (i£U(pi) + \E A (p 2 ) - ±E A (p 3 ) - \E A {p±) - E s ) . (B.12) 

In each of the four collinear limits, two of the E A terms have equal and opposite divergences 
to each other and two of them have equal and opposite divergences to Eg, yielding an 
integrable integrand with a finite result. We have thus proved that Eq. (|B.10[) is finite. 



Although this argument was formulated in terms of the 6-quark axis definition, it applies 
equally well to any infrared-safe definition, like the thrust axis, since they must become 
equal in the triple-collinear limit. 

Since the integrand is everywhere integrable, we can use the same numerical program 
as for the rest of the non-singlet contributions, and obtain the results in Eqs. (|54L |55|). 
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